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Abstract 

We explore kaon condensation using the bound-state approach to the Skyrme model on 
a 3-sphere. The condensation occurs when the energy required to produce a K~ falls below 
the electron fermi level. This happens at the baryon number density on the order of 3-4 
times nuclear density. 
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1 Introduction 

Recently it was discovered |I[ that kaon condensation in nuclear matter of densities on the 
order of a few times that of normal nuclear matter may have an impact on the formation of 
low mass black holes instead of neutron stars for masses on the order of 1.5 solar masses. 

Since the proposal of kaon condensation and its possible consequences in 1986 , various 
attempts have been made to determine the critical density at which kaon condensation sets 
in, by means of a chiral lagrangian p|, ^, ^, ^ and also by means of phenomenological off-shell 
meson-nucleon interactions ^. 

In this paper, we explore kaon condensation within the framework of the bound-state 
approach to the Skyrme model |0, |]. In 1985, Callan and Klebanov proposed the bound- 
state approach to the Skyrme model as a model for exploring the properties of strange 
baryons. To this day, the bound-state approach to the Skyrme model remains one of the 
most successful models for hyperons. At around this time, it was suggested that the SU(2) 
Skyrme model may be used to determine the properties of dense nuclear matter by means of 
putting skyrmions on a cubic lattice [^. It was later determined that these SU(2) skyrmions 
undergo a transition to a smoothed out "dissolved" phase with restored chiral symmetry. 



||10| , |Tl|. The cubic lattice symmetry was also replaced by other symmetries which resulted 
in a lower crystal interaction energy without affecting the qualitative phase structure of the 
crystal fT^.Q 

Unfortunately, lattice calculations can be difficult and require involved numerical work. 
In 1987, Manton proposed an alternative to the lattice. He put a single skyrmion on the 



surface of a 3-sphere of finite radius |jT3[. Various densities of nuclear material can be ex- 
plored simply by varying the radius of the hypersphere. Since there is only one skyrmion, 
the calculation is much more straightforward. The phase transition may be calculated an- 
al5d;ically. While the model seems somewhat contrived, since the nuclear matter density 
is manifested entirely by the curvature of the hypersphere, it nevertheless produces similar 
qualitative and quantitative (to within a few percent) results as the lattice calculations at a 
fraction of the computational cost. Furthermore, a study of two skyrmions on a hypersphere 



IJ] suggests that the simultaneous addition of skymions and the increase in hypersphere 
radius, while keeping density fixed, does not significantly affect the results. Thus, one may 
be able to continuously connect the lattice approach and the hypersphere approach without 
substantial change in the observable quantities. 

In this paper, we explore kaon condensation by adapting the bound-state approach to 



the 3-sphere. This approach has already been carried out in [^ using a different version 



of the Skyrme model, in which vector mesons are included as fundamental fields. It was 



further reasoned in |T^ that the parameter mx, the mass of the kaon in a vacuum, should 
"run" as a function of the density. At the time, this seemed to be the only means in which 
kaon condensation could be derived from the Skyrme model. It was assumed that unless the 



^Unfortunately, kinetic energy effects have not been computed with the Sfcyrme model due to its non- 
renomahzabihty, so such comparison of energies is inherently difhcult. 

1 



energy of the kaon, tu, reaches zero for a finite value of the 3-sphere radius, kaon condensation 
would not occur at all, or at least not for reasonably low densities. Without the running 
kaon mass, u reaches zero only at zero radius. However, with more careful consideration 
we realize that kaon condensation sets in at a point where uj equals the fermi level of the 
electrons [Q, which may be much higher than its rest mass. This reopened the possibility 
of using the original version of the Skyrme model, with no running masses, to calculate 
the onset of kaon condensation. This method has the advantage of requiring fewer input 
parameters, thus pinning down a more precise (albeit model-dependent) transition density. 
It is interesting that the final result lies within the range reported in P| and lies fairly close 
to the values suggested by [|^ . 

This paper is organized as follows: The bound-state approach is summarized in Section ^ 
and adapted for the hypersphere in Section |^. The background field is calculated in Section ^ 
In Section |^, a new discrete symmetry is found and discussed, which leads the way for the 
calculation of the effective kaon mass, a;, as a function of hypersphere radius in Section ^ 
Finally, in Section |^, uj is compared to the electron fermi level and a prediction for the K~ 
condensation transition density is obtained. 

2 The Bound State Approach to the Skyrme Model 

A single hyperon may be modeled by the Skyrme model in fiat space. The version used by 
Callan and Klebanov [^ utilizes a simple chiral lagrangian with a commutator stabilizing 
term 

^ = ^tr(9^f/t9^f/) + ^tr[9^[/t/t,9,[/f/t]2 + ^tr^(t/ + f/^-2) 



■^ tr M^iVf^ + -^tr[M^,M,]2 + 4 tr;W(t/ + t/t - 2) 



2 --L"'"A" "'"J^J ' o --■'•'V'- ' ^ -/) l-'-j 



where U{x,t) G SU(3). Ai is proportional to the quark mass matrix and, if we neglect the 
u and d masses, is given by 



f 







m 



K 



where, in general, we write 3x3 matrices in the partitioned form 



2x2 



1 X 2 



2 X 1 



1 X 1 



M^ is defined as 



M„, := dMUl 



The standard fit to nucleon and delta masses yields f^ = 129 MeV and e = 5.45. The 
empirical value of the kaon mass is mx = 495 MeV. 



In addition to the local terms, we must also include the Wess-Zumino term, which is 
written in non-local form as an integral over a five- dimensional disk with spacetime as its 
boundary: 

^wz = -^^^ J^d''xe^''''''HT{M^M,MaMf,M,), (2) 

where U{x,t) is continuously extended to the disk (with our convention, e'^1234 _ -|^^_ jy jg 
the number of colors. 

In the bound-state approach to the Skyrme model, we use the following anzatz 

Uix,t) = ^U^{x,t)UKix,t)^U^{x,t), (3) 



2^2 
UK{x,t) = expi 



fn \ K{x,ty 



K{x, t) 







where U.,^ G SU(2) and K isa complex spinor. Substituting this anzatz into the Skyrme action 
and expanding to second order in the kaon fields, we obtain the bound-state lagrangian 



C = -lftT{Al) + —tT{[A„A,]')-mlK^K + D,K^D,K 



ftr(4) + ^^ 
+ -^D^K^IA^, A,]D,K - -^D^K^D,KtT{A^A,) 
+ -^T^D.K^D.KtriAl) - -^K^KtT{[A„A,]') + l-K^KtT{Al) 
+ TTTFT^^"^^ triA^ApA,) {D^K^K - K^D^K) , (4) 



where 



The term proportional to A^ is the result of integrating the Wess-Zumino term. The 
lagrangian itself is quadratic in kaon fields and describes non-interacting kaons bound to a 
background SU(2) skyrmion. The background field is realized by the hedgehog anzatz 



Un{x,t) 



piF(r) T-r 



ot 



fl (5) 



with the boundary conditions F{0) = n and F{oo) = 0. We write the eigenmodes for K{x, t) 
in terms of its spin and isospin. In particular, we will be interested in S-wave and P-wave 
kaons, for which the corresponding anzatze are 

Ks{x,t) = ksir)xit), (6) 

Kp{x,t) = ikp{r)f-rx{t). (7) 



In each case, x{t) ^'^'^ the dynamical variables arranged in a complex spinor. 

The field F{r) is determined classically using the non-kaon part of the lagrangian. This 
baryon field is treated as a background field with which the kaons interact. The energy 
eigenstates form a Fock space for each kaon mode, where the energy required to create a 
kaon in a particular mode is given by its classical eigenfrequency (i.e. substitute x = —i^X 
and minimize the action with respect to u and k{r)). In fiat space, the lowest energy modes 
for P-wave and S-wave are given respectively by up = 153 MeV and cus = 368 MeV. 

3 The Skyrme Model on a Hypersphere 



In |T^, Manton suggested that a baryon crystal may be approximated by putting a single 
Skyrmion on a hypersphere of finite radius. This has the effect of allowing the baryon to 
interact with itself. We will continue with this idea by considering the bound state between 
a kaon and a baryon on the hypersphere. 

A hypersphere of radius a is described by three angular coordinates (p, 0,(f)) {0 < p,9 < it 
and < < 27r) with the following spatial metric 

ds"^ = a^{dp^ + sin^ pidO"^ + sin^ Odcj)'^)). 

It is often convenient to replace the p coordinate with x = cot p, — oo < x < oo. This spreads 
out the coordinate system around the poles, thereby eliminating coordinate singularities in 
the lagrangian. 

Topologically, the hypersphere is equivalent to R^ with all the points at infinity mapped 
to a single point. However, the metric is different, with points near infinity (r ^ a) much 
closer together on the hypersphere than the corresponding points in fiat space. For r <^ a, 
this difference is negligible. Since the baryon size in flat space is on the order l/e/jr, it follows 
that for a ^ l/ef^^, the results on the hypersphere will not be much different from the results 
in flat space. However, for a on the order of l/e/jr, the skyrmion begins to interact with 
itself on the hypersphere, and this can have a significant effect on the results. 

We choose to identify the origin r = in R with the south pole p = vr (x = — oo) and the 
point at infinity with the north pole p = {x = oo). For r <ti a, S^ and R^ are equivalent, 
both in topology and metric, with corresponding points 

r ~ a{7T — p) ^ —a/x. 

There are a few technical remarks in adapting the Skyrme model for the hypersphere. 
The profile functions F and k on the hypersphere will depend on p (x), which replaces r as 
the radial coordinate. The boundary conditions for F are as follows 

at the origin at infinity 



r 





OO 


P 


IT 





X 


— OO 


oo 


Value of F 


TT 






In choosing to identify the origin with p = it, we reverse the p direction compared to the 
f direction in flat space. It is therefore appropriate to reverse the sign of £t^'^°'f^^ and hence 
the Wess-Zumino term.Q In flat space, the matrix r ■ f depends on the angular coordinates 
6 and 0. On the hypersphere, we identify r ■ f with the same function of 6 and 0, except 
these coordinates are now the hypersphere coordinates 6 and 0. 

Using the anzatz in (||) adapted for the hypersphere, we may compute A^ and V^. In 
terms of the normalized coordinate basis for the dual tangent space 



u" = dc^^\gcc\, c = t,p,9, 
we find V^ and A^ to given by 



1 ^,. ^ ^ _ „ sinF 
-F'pzr-fa;^ + ^- 
2a ^ 2r 



A, = _F'(p)zr-f.;^ + ^^zr-(^c^J + 0.;5 



where r = rphys = a sin p. Using this result and (|^), we obtain the (spatially integrated) 
lagrangian for a P-wave kaon 

tt/^ /■°° , f /sin^F 1 



e J-oo 

1 

+-- 
a 



, f /sin^ F I , ^, , ..2\ 
2 sin^ F + 4(1 + x^) sin^ F{F' (a;))' 



efi, J-oo [ \ a 

24 
+ — (1 + a;^) (1 + cos F) sin F F' ix) k ix) k' (x) 
a 

\(l+x^) ^ l + x-^^J 



AeujN 



+ 4^ sin^ F F' (x) (A; {x)f + (k {x)f (a ( - (F' (x))' 



(1 + cosF) sin^F\ 8/1, n, , ^, o/^// nn: 

+ 2 ^ ^ , ' - 2—^) + -(-1 + x^ 1 + cosFf F' X y 
1 + x2 1 + x2/ a V4' 

- 2(1 + x^) sin^ F (F' (x)) V (1 + cos F)^ sin^ F - sin^ F 
+ ^ 



eV/ (l + a;2) 

where a = aef-,, and we use the substitution x = cotp mentioned earher. One could use 
(P) instead of (|^ to obtain the lagrangian for S-wave kaons. However, as we will see later, 

^If we choose instead to identify the origin with p = 0, then F(0) = tt and F{Tr) — and the Wess-Zumino 
term stays the same. The resuhs are equivalent. 



it is possible to obtain this lagrangian from the P-wave lagrangian by substituting F{x) i-^ 
F{x) + TT. The ramifications of this result will be explored later. 

The remainder of the paper will be spent solving the stationary conditions for F{x) and 
k{x), and determining u as a function of a. 

4 Finding the Background SU(2) Baryon Field 

The SU(2) baryon field is determined by the profile function F{x), which is obtained by 
solving the stationary conditions for F using the non-kaon part of the lagrangian in (§). 
Much of this work has been carried out by Manton [0 using different methods. 
The differential equation for F{x) is given by 

1 „„ , , cosFsinF 4 /, ,, n „ „„ , , 
= -F"{x) -—— + —({l + x^)sm^FF"{x) 

2 1 + x"^ a"^ ^ 

+ (1 + x^) cos F sin F {F' {x)f + 2x sin^ F F' {x) 

-cosFsin^F). (9) 

The boundary conditions for F are given by 

F(-oo) = 7r, F(cx))=0. 

The identity profile function F{x) = arccot(x) {F{p) = p) is a solution of (H) for every 
a. It corresponds to a baryon uniformly distributed on the hypersphere. For a less than the 
critical value acr = 2\/2 it is the only solution and it is stable. 

As a crosses the critical value, the identity solution becomes unstable in favor of an 
asymmetric solution, where the baryon is located at one of the poles. There are actually two 
equivalent solutions which correspond to each other under the transformation 

F{x) H^7r-F(-x), 

which leaves the differential equation invariant. We always choose the solution with the 
baryon localized around x = — oo (where F = 7r).0 

To find the localized solution, we numerically solve the differential equation using the 
initial conditions provided by the asymptotic solution 

F(x) ^ vr + A/x + ... X < 

and use trial and error with the value of A until we find a solution which approaches zero as 
X -^ oo. For a ^ ^/ef-n (i-e. a ^ 1), we expect F{x) to approach the solution corresponding 
to fiat space. In particular, F takes its middle value (7r/2) at around r = asinp ~ I/c/tt 
which corresponds to x ~ —a. Thus A ^ a for large a. 



•^Localizing the baryon around the F = n pole, as opposed to the F — pole, results in a slightly lower 
energy for the physically relevant case to-tt 7^ 0. 
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5 A Discrete Symmetry 

In section ^ we found that the non-kaon part of the lagrangian is symmetric under the 
transformation 

F{x)^n-F{-x). (10) 

This transformation maps the baryon from one pole to the other. A simple coordinate 
transformation corresponding to point inversion, which also maps the baryon from one pole 
to another, is given by 



t-^ + vr, 

e ^ TT-e, 

X H-> —X, 



F{x) 



-F(-x), 



;ii) 



and leaves the original Skyrme lagrangian invariant .0 (|T0]) differs from the coordinate inver- 
sion of ([llD by the transformation 



Fix) 



F{x) 



TT. 



(12) 



If one examines closely, one realizes that the kaon part of the P-wave lagrangian in (§) 
is not symmetric under ( p!oD or ([T2|) , regardless of the substitution for k{x), although it is 
symmetric under (|ll]). However, (p!2D (and hence ([lOD) corresponds to a valid symmetry of 
the bound-state lagrangian. It is no accident that the non-kaon part is symmetric under 
(p!0|). To determine what this symmetry transformation does to the kaon field, we need to 
examine the effect of (|12[) on the Skyrme lagrangian. 

The transformation in (|l^) corresponds to 



U^ 



ot 







f/.. 



The bound-state approach requires the value of ^/U^. Technically, the square root is am- 
biguous because [/^ is not positive-definite hermetian. We define the square root to be 



■ F(x) 



ot 







Hence (|1^) corresponds to 



IT ■ r 



Ot 



1 



U^ = \U^ 



IT ■ r 



Ot 



*The boundary conditions on F{x) are different. 



If we arrange for the kaon field to transform so that 



Uk ^ W^UkW, W : 



then U{x, t) itself will transform as 



tT ■ r 



ot 



(13) 



U 



Ot 



1 



u. 



The Skyrme lagrangian (|T]) is invariant under ([Mf) . The transformation in (|n 
to 

K I— i> — ir ■ rK. 



(14) 



corresponds 



(15) 



If we compare the anzatze in (^ and (^ we find that, aside from a trivial change in the sign 
of k{x), (|T5|) corresponds to the exchange in the roles of S-wave and P-wave. The crucial 
conclusion from the discussion above is that the simultaneous transformations ( PUJ ) and (|T3p 
leave the action invariant .[| This simultaneous transformation does the following: 

• It moves the baryon from one pole to the other. 

• It does not affect the boundary conditions of F{x). Therefore, the value of F at the 
pole which the baryon is nearest to changes from vr to or vice-versa. 

• It transforms the kaon from S-wave to P-wave, but does not otherwise move the kaon 
{k{x) does not change, except for an irrelevant change in sign). 

For a > Oct the solution for F{x) does not respect (p^. In this phase, the kaon field is not 
expected to respect (pISj) either. Thus, there is no symmetry between the S-wave and P-wave 
kaons. However, for a < «„, the solution F{x) = arccot(a;) (-F(p) = p) does respect ([T0|). 
In this phase, there is no distinction between S-wave and P-wave kaons. In the next section, 
we see the value of u for S-wave and P-wave approach each other as the critical value of a 
is approached. For a < Ocr, there is only one value of u. 



6 Solving for the Kaon Field 

The next goal is to solve the kaon part of the lagrangian in (^ for uj and the kaon field 
as a function of a. In the next section, we will compare the electron fermi level to u and 
determine the value of a where kaon condensation sets in. 

In the non-uniform phase a > 2-\/2, we substitute the numerical solution to F(x) into 
the Euler-Lagrange equations corresponding to the lagrangian and numerically solve for k{x) 
and uj. We do this separately for P-wave and S-wave. 



^This symmetry is broken by introducing a pion mass, which we neglect in this paper. 



For P-wave kaons, we use the lagrangian given in (|]). The initial conditions are provided 
by the asymptotic solution 

k{x) ^ fco + • • • X <C 0. 

Since the differential equation is linear in k{x), we may fix the value of kg = 1. The value 
of uj is adjusted until we find a solution which approaches zero as x ^ oo. A typical 
(unnormalized) solution for k{x) is given in Fig. |I|. 

For S-wave kaons, we use the S-wave lagrangian obtained from (H) by the substitution 
F{x) I— > F{x) + TT. This transformation differs from ([T0|) by coordinate inversion, and so 
the kaon switches poles instead of the baryon. The initial conditions are provided by the 
asymptotic solution 

k{x) ^ -A/x + ... X < 0. 

We fix y4 = 1 and adjust the value of u until we find a solution which approaches a constant 
as X — s> oo. A typical (unnormalized) solution for k{x) is given in Fig. |^. 

The solution for u as a function of a for both S-wave and P-wave is given in Fig. ^. 

In the uniform phase a < 2-\/2, we substitute the analytic solution F(x) = arccotx 
{F{p) = p) into @. As explained in the previous section, there is no distinction between 
S-wave and P-wave kaons in the uniform phase. If we substitute the same analytic form for 
F{x) into the lagrangian appropriate for S-wave, we obtain the same differential equation, 
except for the coordinate inversion x ^-^ —x. 

The resulting differential equation is 

k{x) 
where 



k"{x) + ,^ ; J (- - xVl + x^ - x' + A) = 0, (16) 



_ 3 + a^ / uja\ e^N ua _ 1 (mK\ 4 . . 

~2 + a^\eU) 7c^{2 + a^)eU 2 + a'[efJ "' ^ ^ 

The dependence of cj on a is entirely contained in the expression for A. In the intro- 
duction, we claimed that u; — > as a — > 0. From (|1^, it is clear that this is only possible 
if there exists a solution to (|16]) for A = 0. This is not a trivial prediction and there is no 
reason to expect it. However, if we rewrite (0) in terms of p instead of x, we do in fact find 
a solution for A = 0. It isQ 

k{p) = sin -p. (18) 

It is this miracle which ultimately is responsible for the fact that the Skyrme model yields 
kaon condensation at reasonable densities, if at all. Otherwise, u ~ A/a, which is the same 
a-dependence as the electron fermi energy. It would not be clear for what values of a, if any, 
that u < Ep. 



^This solution was also found in a different version of the Skyrme model pi 



Solving (|T7|) with A = for cj as a function of a, we findQ 



^a _ V ^S^)' + (^)'"'(3 + a^)-g ^ ^^g^ 



e/,r 3 + a^ 

After plugging in numbers, we get 



'I + .02433a4(3 + ^2) _ i 
u = (3174 MeV)^^ j—^ '- . (20) 

At the critical point, we find uo = 332 MeV, which is consistent with the numerical solutions 
for a > 2-\/2- For a — i> 0, we find 

uj ^ (38.6 MeV)a='(l + order " ^^J" " ^ ). (21) 

It should be noted that k{p) given in (0) satisfies 

k{p) = sin -F{p). 

This is the form for k{p) which corresponds to a small rigid rotation into the strange directions 
0. Thus the rigid rotator approximation to the bound-state approach given in |]16| becomes 



exact in the uniform phase. Indeed, the formula (|T9]) for ui can be expressed in the form 



obtained in [16| using the rigid rotator approximation: 






8$ VV Mo' J 4vr$ 

The appropriate values for $ and F are 

7r2a(3 + a2) 



TT^a^ 



^^ Jtt ^ Jtt 



$ 

In light of the fact that k{p) as given in ( |I8| ) is independent of mx, it is no surprise that 
the form of k{p) matches the rigid rotator form. For rriK = 0, we expect the SU(3) flavor 
symmetry to be restored. 

The only remaining question is whether or not A = actually yields the lowest value of 
uj for all a below the critical value. We don't expect A 7^ to yield lower energy solutions 
for any a < 2-\/2 because 

• uj{aci) is predicted correctly by A = 0. Also, the A = solution for k{x) matches the 
numerical solutions we obtain for k(x) for a > 2v2 for both S-wave (inverted) and 
P-wave. 



^The other solution to the quadratic equation (corresponding to a minus sign in front of the square root 
sign) corresponds to lu < which represents the energy of a K'^ in nuclear matter. As expected, |w| diverges 
as N/a as a ^- due to the Wess-Zumino term. 
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• k{p) given by the A = solution contains no nodes, which is what we expect for the 
ground state solution. 

• k{p) given by the A = solution matches the rigid rotator solution, a situation we 
expect for rriK = 0. 

• A = must yield the lowest u for a sufficiently close to zero. All other values of A 
result in a divergent uj as a ^ 0. 

• A sudden jump from A = to a different value of A would involve a discontinuous jump 
in u and in k{p), neither of which is expected. 

7 Determining the Onset of Kaon Condensation 

To study the onset of kaon condensation, we need to understand what is going on in the 
neutron star. In our model, the neutron star consists of a Skyrme lattice of protons and 
neutrons, as well as enough electrons to balance the charge of the protons. For simplicity, we 
shall treat these particles as free particles. We assume that the star is in quasi-equilibrium. 
In particular, the following reaction 

is in equilibrium. This together with conservation of baryon number and charge allows us 
to solve three equations 



tp - 


= 4 + 6^, 


PB = 


= Pp + Pn, 


Pp -- 


= Pe, 



for the three variables p„, Pp, Pe, in terms of a. The baryon number density is given by the 
inverse volume of the hypersphere 

PB = {2n^a^)-\ 
The fermi energy €p for a spin 1/2 fermion is given in terms of its number density by 



£f = ym^^+kj, 



^m) + (3vrV/)2/='. (22) 



This will allow us to calculate ey as a function of a. Before doing so, we wish to discuss 
kaon condensation. In free space, the conversion of electrons to kaons via the reaction 

11 



is unthinkable, because the mass of the kaon is much larger than that of the electron. How- 
ever, from the previous section we see that the effective mass of the kaon, cj, approaches 
zero for large baryon densities. Furthermore, for large baryon densities, the fermi level of the 
electrons (as well as the protons and neutrons) increases. It follows that there exists a value 
of a for which e'^p = uj. Beyond that density, it becomes favorable for electrons to become 
kaons. Since kaons are bosons, there is no fermi level for the kaons. Every kaon requires 
the same energy lo to produce. As explained in |I|, kaon condensation on the order of a 
few times normal nuclear matter density may have an impact on the formation of low mass 
black holes, by reducing the resistance of the neutron star to collapse. We now compute the 
density of kaon condensation as predicted by our model. 



The formula in (|22D for the fermi energy is fully relativistic. We shall make the approx- 
imation that the fermi-level electrons are ultra-relativistic and the fermi-level nucleons are 
non-relativistic. Thus the fermi energies are given by 

Conservation of charge and baryon number allows us to express p„, pp and pe in terms of 
one parameter, 7 = Ppj Pb- 

Pp = Pe = IPB, p„ = (1 - j)pb- 
In terms of this parameter, the fermi levels may be reexpressed 



Sp - 


= {'iTr^pBY^'Pil-lf' + mj, 


4 = 


- (3vrVi.)i/3/3f/3 + m^, 


^^ = 


= {^7:^pBf'W\ 


/3 := 


{Sn^PBy/' 



iN 



Plugging in numbers yields 



{Stt^PbY^^ = (805/a)MeV, 
p = 0.429/a. 

The equilibrium equation then becomes 

/3((1- 7)^/3- 7^/3) =y/^ (23) 

which we can solve for 7 in terms of a. 

To get an order of magnitude estimate, we calculate the fermi energies at a = 2-\/2- This 
corresponds to f3 = 0.152. Solving (El) iteratively for (3 small, we find 



7 = /^•^ - 3/3^ + . . . f« 3.27 X 10 
12 



-3 



from which we compute the fermi levels 



ep{(^cr) = 


= 42.2 MeV, 


^^("cr) - rriN = 


= 0.95 MeV, 


Spioicr) - rUN = 


= 43.2 MeV. 



The relativistic approximations made earlier are justified by these results. Note that ti;(Q;cr) = 
332 MeV is much larger than Sp at the phase transition, so that kaon condensation must take 
place for a < 1\pl. However, oj decreases quickly from its critical value and e%^ continues 
to rise.0 For a = 1, cu ~ 40 MeV, which is definitely lower than the fermi level at a = 1. 
Therefore, kaon condensation occurs somewhere in the regime 1 < a < 2-\/2. 

To find the transition exactly, we note that (^) can be explicitly solved for a in terms 
of 7: 

fl _ ^)2/3 _ 2/3 

Hence both sides of the equation e^^ = lo can be written explicitly in terms of 7. Numerical 
solution yields 7 = 0.016.0 The following values are found at the onset of kaon condensation 



7 


= 


0.016, 


OL 


= 


1.58, 


a 


= 


0.44 fm. 


UJ = Ep 


= 


129 MeV, 


£p -rriN 


= 


8.8 MeV, 


Ep — rriN 


= 


138 MeV, 


PB 


= 


0.595 fm"^ 



The density for kaon condensation is on the order of a few (3-4) times the density of normal 
nuclear matter pnuci = 0.16 fm~ . 

8 Conclusions 

We have predicted the onset of kaon condensation to be pB = 0.595 fm~^ = 3.7pnuci- We offer 
a brief comparison to the literature in Table [^. Beyond the critical density, an increasing 
fraction of neutrons is converted to protons and K~ , which softens the star's resistance to 
compression. As explained in ||l|], a critical density on the order of a few times the nuclear 

^In the uniform phase, the nuclcons lose their identity, so the calculation of the fermi level might seem 
suspicious in this regime. Nevertheless, it is reasonable to expect the fermi level to continue to rise and 
for (p2) to be a reasonable approximation. 

^This means protons make up 1.6 percent of the nucleons in the neutron star. Most of the protons have 
already recombined with electrons to form neutrons. 
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Table 1: Values of the baryon transition density pb 

Transition density Relative density 

Source Pb Ps/Pnud 

Bound-state approach 0.595 fm^'^ 3.7 

on a hyper sphere 

Mean field approach 0.3-0.6 fni~ 2-4 

with chiral lagrangian P 

Vector-meson Skyrme 0.15-0.23 fm^ 0.9-1.4 

model on hypersphere ||T5| 



density is sufficient to reduce the upper mass limit beyond which the collapse of a star to a 
black hole is inevitable. Our results fall within that range. 

The key ingredient which allows the bound-state approach to predict a reasonable tran- 
sition density is using the electron fermi level. Without the fermi level, kaon condensation 
would not set in until infinite density. 

Now that it has been shown that the Skyrme model can produce reasonable results for 
K'^ condensation, several refinements of the calculation are possible. Since the rigid rotator 
approximation to the bound-state approach becomes exact for a < 2\/2, we may calculate 



the complete 1/A^ correction to the K~ energy using the formulae generated in [T^. A 
more ambitious project is to replace the hypersphere with a real lattice. In addition, a more 
sophisticated model of the neutron star may be used to determine the electron fermi energy 
more precisely. None of these refinements is expected to drastically change the results. 

Acknowledgements 

I would like to thank I. R. Klebanov, my thesis advisor, for recommending this problem to 
me and for sharing his ideas. This work was supported in part by DOE grant DE-FG02- 
91ER40671, the NSF Presidential Young Investigator Award PIIY-9157482, and James S. 
McDonnell Foundation grant No. 91-48. 

References 

[1] G. E. Brown and H. A. Bethe, Astrophys. Jour. 423 (1994), 659. 

[2] D. B. Kaplan and A. E. Nelson, Phys. Lett. B175 (1986), 57. 

[3] G. E. Brown, K. Kubodera, M. Rho, Phys. Lett. B192 (1987), 273; H. D. Politzer, M. B. 
Wise, Phys. Lett. B273 (1991), 156; G. E. Brown, K. Kubodera, M. Rho, V. Thorsson, 



14 



Phys. Lett. B291 (1992), 355; V. Thorsson, M. Prakash, J. Lattimer, Nordita-93/29 N, 
SUNY-NTG-92-33. 

[4] G. E. Brown, C.-H. Lee, M. Rho, V. Thorsson, Nucl. Phys. A567 (1994), 937; C.-H. 
Lee, H. Jung, D.-P. Min, M. Rho, Phy.s. Lett. B326 (1994), 14; C.-H. Lee, G. E. Brown, 
M. Rho, SNUTP-94-28, |hep-ph/9403339i 

[5] C.-H. Lee, G. E. Brown, D.-P. Min, M. Rho, Preprint |hep-ph/94063TIl (June 1994). 

[6] H. Yabu, S. Nakamura, F. Myhrer, K. Kubodera, Phys. Lett. B315 (1993), 17; M. Lutz, 
A. Steiner, W. Weise, "Kaons in dense matter", to be pubhshed. 

[7] C. Callan and L Klebanov, Nucl. Phys. B262 (1985), 365; C. Callan, K. Hornbostel 
and L Klebanov, Phys. Lett. B202 (1988), 269; L Klebanov, Strangeness in the Skyrme 
Model, in Hadrons and Hadronic Matter, Plenum (1990). 

[8] N. Scoccola, H. Nadeau, M. Nowak and M. Rho, Phys. Lett. B201 (1988), 425; J. P. 
Blaizot, M. Rho and N. Scoccola, Phys. Lett. B209 (1988), 27; N. N. Scoccola, D. P. 
Min, H. Nadeau and M. Rho, Nucl. Phys. A505 (1989), 497; M. Rho, D. O. Riska and 
N. N. Scoccola, Z. Phys. A341 (1992), 343. 

[9] I. Klebanov, Nucl. Phys. B262 (1985), 133. 

[10] E. Wust, G. E. Brown, A. D. Jackson, Nucl. Phys. A468 (1987), 450. 

[11] A. D. Jackson, J. J. M. Verbaarschot, Nucl. Phys. A484 (1988), 419. 

[12] M. Kugler, S. Schtrikman, Phys. Lett. B208 (1988), 491. 

[13] N. S. Manton, Comm. Math. Phys. Ill (1987), 469. 

[14] A. D. Jackson, A. Wirzba, L. Castillejo, Phys. Lett. B198 (1987), 315. 

[15] H. Forkel, A. D. Jackson, M. Rho, N. N. Scoccola, Nucl. Phys. A509 (1990), 673. 

[16] D. Kaplan and I. Klebanov, Nucl. Phys. B335 (1990), 45. 

[17] K. Westerberg and L Klebanov, Phys. Rev. D50 (1994), 5834. 

Figures 

Fig. |l] P-wave kaon wavefunction (unnormalized) for a = 5. 
Fig. |2| S-wave kaon wavefunction (unnormalized) for a = 5. 
Fig. ^ Kaon energy (MeV) vs. a. 
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Figure 1: P-wave kaon wavefunction for a = 5 



Figure 2: S-wave kaon wavefunction for a = 5 



Figure 3: kaon energy (MeV) vs. a 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9411430vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9411430vl 



This figure "figl-3.png" is available in "png" format from: 
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